A mechanism for air induced fiuidization in vibrated granular beds 
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We present a new mechanism for the fiuidization of a vibrated granular bed in the presence of 
interstitial air. We show that the air flow induced across the bed, as the gap in the bottom of the 
cell evolves, can fluidize the bed in a similar way as in gas-fluidized static beds. We use the model 
of Kroll to quantify the relevant variables suggested by the mechanism proposed. The relevance of 
the above fiuidization for segregation phenomena is discussed. 

PACS numbers: 45.70.-n, 47.55. Kf, 45.70.Mg, 64.70.-p 



A forced collection of grains can exhibit very rich and 
striking behavior that sometimes resembles the proper- 
ties of fluids and solids 0, 0- Fluidlike behavior can 
be induced, among other means, by shaking @> or 
by injecting gas through the grains 0, IE U ■ Due to 
its particular transport properties gas-fluidized granular 
beds have found many aplications in industry |4| . Vibra- 
tion commonly induces undesirable segregation in many 
industrial processes which involve the handling of large 
collections of grains 0, • From a fundamental point of 
view, forced granular beds offers a plethora of interest- 
ing collective behavior. Thus a basic understanding of 
phase changes as a response of these dissipative systems 
to a change in forcing conditions, grains properties and 
of any relevant parameter is of practical and fundamental 
interest. 

The influence of interstitial air in phenomena related 
to vibrated granular media have been studied in different 
contexts IE 0, E3, E3, E3, E3l Under certain condi- 
tions, interstitial air has proved to be a crucial ingredient 
in segregation phenomena 0, 0, 0, 0, Q , and in this 
line some efforts to include air in molecular dynamics 
simulations have appeared recently |l3j . However, a co- 
herent description that would allow us to understand the 
role played by the air in vibrated granular systems is still 
lacking. 

In this article we present a mechanism by which a vi- 
brated bed of grains can become fluidized in the presence 
of interstitial air. The key ingredient is not just the pres- 
ence of air but the way it flows through the vibrated bed. 
We discuss the consequences of this kind of fiuidization 
in terms of a previous ly p roposed quantitative model [TEj 
for reverse buoyancy [l6| . 

We assume that the vibrated grain pack acts like 
a porous piston so the acceleration of the container 
induces a pressure difference between its ends as the 
pack moves relative to the bottom of the container, 
which in turn produces an air flow across the bed. We 
use Kroll's model to describe the evolution of the gap 
under the grain pack, within one period of vibration. 
This model takes into account a force that arises as a 
consequence of the pressure difference across the bed 
that develops as the gap in the bottom evolves. We 
propose that the flow of air induced by this pressure 
difference can fluidize the medium in much the same 



way as it occurs in gas-fluidized static beds, where gas is 
injected against gravity. The model of Kroll will allow 
us to present explicit relations between the flow of air 
through the bed, the pressure difference between its ends 
and the gap formed at the bottom of the bed. From 
these relations we develop a quantitative description of 
the proposed mechanism of fiuidization. 

Kroll's Model 

In the Kroll's model the granular bed its treated 
as a porous piston. When the effective gravity starts to 
point upward the piston takes off. According to Kroll's 
model, in the nonincrtial reference frame placed on the 
container, only two forces are taken into account during 
the flight: the effective weight of the granular bed and the 
force due to the pressure difference across the piston. The 
following evolution equation for the gap at the bottom of 
the bed is obtained (see figure ^): 
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where s is the size of the gap, p a is the atmospheric pres- 
sure, pi is the time dependent pressure in the gap, A is 
the transversal area of the container, m is the mass of 
the piston, g is the acceleration of gravity and w is the 
position as a function of time of the container. 

We can use an equation of state to relate the gap, and 
hence the volume of the air in the gap, with the pressure 
Pi. The air is assumed to behave as an incompressible 
fluid ■ This results in the additional equation 
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where G is the mass of air in the gap and p a is the density 
of air. It is assumed that the flux of air throught the bed 
is given by Darcy's law 0EJ which gives us the following 
equation 
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where k is the permeability of the piston, h it's height 
and n is the viscosity of air. Equations iJTJ-JSJ can be 
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of s and is given by 



w=a sin[oo(/+? )] 



FIG. 1: The Kroll's picture: the granular bed is treated 
as a porous piston. The container is subjected to sinusoidal 
vibrations of frequency u and amplitude a. At the top of 
the bed we have atmospheric pressure p a and Pi is the time 
dependent pressure in the gap. 



integrated for the unknowns s(t), Pi(t) and G(t). We 
assume that the container is subjected to sinusoidal os- 
cillations of amplitud a and angular frequency ui, so we 
have that w = a sin[w(t + to)], where sin(w^o) = 1/r and 
r = aw 2 /g. At t = the medium takes off. 

If we introduce the adimcnsional quantities s' = s/a 
and t' = tot, we find the following evolution equation for 
the adimcnsional gap s': 



s' + —s' = i[r s in(t' + 4) 

T K I 



1]. 



(4) 



where s' = ds' jdt 1 , t' = coto and the adimensional num- 
ber tk is given by 



tk 



ujkp Tl 



(5) 



We see that in this model the effect of air in the piston 
movement is a usual frictional factor, with a friction coef- 
ficient given by 1/tr-. We expect then that if tk is small 
then the gap will be small. The adimcnsional number tk 
depends on the grains diameter through the permeabil- 
ity k of the medium |17| . An interesting prediction that 
arises is that the gap's dynamics depends only on the 
adimensional numbers T and tk- Integrating equation 
(0J with initial conditions s'(0) = s'(0) = we obtain 
the evolution of the adimensional gap s'(t'): 

s' = 2£ [C-f-T K (1 + AT sin S)e- t,/TK 

-\Tcos(t' -8)], (6) 

where C = tk + AT (tk sinJ + cos 6), A = (l + t k ) , 
and S — arctan tk — arcsinl/T. Equation lO is valid 
while s' > 0. 

The gap is small 0] when a large pressure difference 
between the ends of the bed is induced. From equations 
@ and we see that in this model the pressure differ- 
ence across the bed is proportional to the first derivative 



Pa ~ Pi 
ghp m 



(7) 



From Darcy's law, the flow of air through the grains 
is proportional to this pressure difference. 

Gas-Fluidized beds and cyclic fluidization in vibrated 
granular beds 

In the technique of fluidizing beds by injecting gas 
against gravity a fluidized state is achieved when the 
force due to the pressure difference across the bed is high 
enough to cancel the weight of the bed 0. 0- Q An 
important control parameter in this technique is the ve- 
locity U at which gas is injected upward at the bottom of 
the static bed. As the gas flows up through the granular 
medium it induces a pressure difference that results in an 
effective force on the bed that point upward. This force 
rises as U is increased from zero and when it equals the 
weight of the bed fluidlike behavior starts to be observed 
14 |5l . This fluidization condition can be written as 



APA > mg, 



(8) 



where AP is the pressure difference between the ends of 
the bed. 

In the Kroll's model the pressure difference and the 
effective gravitational field are time dependent. This can 
be seen if we rewrite the evolution equation Q for the 
gap as 



m ^2 = m 9ef(t) ~ [Pa ~ Pi(t)]A, 



(9) 



where g e f(t) — g{Tsm[co(t + to)] — 1} is the effective 
gravitational field. Before taking off, when g e j points 
downward and the bed is in contact with the bottom of 
the container, there is no significant air flow, the bed is 
being pushed downward against the bottom of the con- 
tainer and there is no relative movement between the 
grains; under these conditions we consider the bed as be- 
ing in a solid state. When g e f starts to point upward the 
bed takes off, the size of the gap and the air flow through 
the grains starts to increase. Following condition (JSJ, for 
this case we can write the fluidization condition as follows 
(9ef + 0): 



[ Pa -pj(t)]A 

mg ef (t) 



> 1. 



(10) 



Combining equations @ and (|10|) , the fluidization con- 
dition H1U|) is initially satisfied at the inflection point of 
the curve s' vs t' shown in figure |21 This result is in- 
dependent of the validity of equations and (J3J and, 
consequently, of the assumptions made in deriving them. 
If we assume that equations © and J^J holds, the be- 
ginning of the fluidized state occurs at the adimensional 
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time t'f (measured from the time when g e f starts to point 
upwards) which is given by 
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Within the Kroll's assumptions, at t'f the pressure dif- 
ference across the bed is maximum. From equations 
@ or and the discussion above, it is expected that 
the mechanism proposed will be effective if tk is small. 
Equation gives, to our knowledge, a new result. 

With this equation we can estimate when the new mech- 
anism of fluidization is initiated. This can be experimen- 
tally verified by measuring what happens to the granular 
medium at the inflection point of the curve that repre- 
sents the evolution of the gap. This, for example, predicts 
for the case of reverse buoyancy, a change in the motion 
of an intruder at time t'^ . 

In figure [21 we show the evolution in time of the gap 
(from eq. ©), the pressure difference across the bed 
(from eq. 10), the effective gravitational field acting on 
the bed and the position of the container for the experi- 
mental conditions of reference [l^. Within the fraction 
of the cycle that runs from t'j until the bed hits the 
bottom of the cell at the end of its flight, we can idcntific 
four zones depending on the direction and intensity of the 
effective gravity g e f and of the air flow. The fluidization 
condition (|10|l is satisfied in zones I and IV (see figure). 
If the flight time of the bed is small compared with the 
period of vibration, then we have a granular bed that 
alternates cyclically between a solid and a fluidized state. 

Cyclic Fluidization and Reverse Buoyancy 

We are now going to explore the consecuences of the 
above mechanism of fluidization in the context of the 
phenomenon of reverse buoyancy. Under certain condi- 
tions large heavy objects immersed in a vibrated granular 
medium rise and similar light ones sink to the bottom. 
This was called reverse buoyancy |l6| . In reference |l5( it 
was proposed that this phenomenon could be explained 
by assuming that the granular bed behaves as a fluid in 
only a fraction of the cycle of vibration. The necessary 
feature is that during most of the time in which the bed 
is in a fluidized state the effective gravity points upward. 

In the fluidization condition explained in the previ- 
ous section, the air flow points downward just after the 
bed takes off, while the effective gravity points upward 
[l9j . If a buoyancy force arises in the above fluidized 
state, large heavy objects immersed in it may sink up- 
ward while similar light ones may float downward. This 
case is portrayed in zone I shown in figure In zone 
IV the fluidization condition IjlUI) is also satisfied, but 
now g e f is pointing downwards and we expect to observe 
regular buoyancy in this zone. Further work is needed 
to anticipate how the bed would behave in zones II and 
III, but it is expected that relative movement between 
the grains is posible in these zones. In particular, we can 
ask ourselves: is there a buoyancy force in zones II and 
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FIG. 2: (Color online) Within a period of vibration, we show 
in this figure the evolution in time of the following adimen- 
sional variables: the gap s' (from eq. ©), the pressure dif- 
ference across the bed (pi — p a )/gh.p m (from eq. 0), the 
effective gravitational field acting on the bed <? e / /<? and the 
position of the container w/a for the experimental conditions 
of reference [T^l . The begining of the fluidized state occurs 
at the inflection point of the curve s' vs t' . It is shown in 
which fraction of the cycle the bed is considered to be in a 
solid state. Within the fraction of the cycle that runs from 
t'f until the bed hits the bottom of the cell at the end of its 
flight, we identify four zones depending on the direction and 
intensity of the effective gravitational field and of the air flow. 
The direction of these quantities is shown at the top of the 
figure. 



III?. In this line, valuable information can be obtained 
by resolving the motion of an intruder immersed in the 
bed within a period of vibration and comparing it with 
the evolution of the measured pressure difference across 
the bed. Other mechanisms can play a role in these zones 

The time that the bed would spend in each of the zones 
shown in figure[2]will depend on experimental conditions. 
Within the model of Kroll, it would depend on T and tk ■ 
If reverse buoyancy is observed, then it could be a sign 
that zone I prevails. In figure we show the time that 
the bed would spend in zones I (Tj) and IV (Tjv) as 
predicted by the model of Kroll. It can be seen that if 
tk is small then we have an scenario that favors reverse 
buoyancy. The model predicts that there is a particular 
value tk,c(X) for which Tj — Tjv- For tk — tk,c the 
average motion of an intruder immersed in the bed would 
depend on how the bed behaves in zones II and III. 

In reference |12| reverse buoyancy was supressed by 
evacuating the air of the system. If we quench the flow 
of air through the bed for the experimental conditions of 
reference [l5| by keeping the top of the container open 
and introducing a filter permeable to air in it's bottom, 
as was made in [fl| , we do not observe reverse buoyancy: 
light and heavy intruders rise. In both procedures the 
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FIG. 3: (Color online) Time that the bed will spend in zone 
I (Ti) and in zone IV (Try) of figure |2] as a function of tk 
calculated from the model of Kroll, for V — 2 and V = 5. 
Ti and Try are reported as its porcentual fraction with re- 
spect to Ti + Tii + Tin + Tiv, that is the time in which an 
intruder immersed in the bed can move with respect to the 
small grains. For tk > 3 the flight time of the bed is close to 
the period of vibration. 



proposed fiuidization mechanism cannot take place. 

In this article we have presented a mechanism by which 
a vibrated granular bed can be fluidized in the presence 
of interstitial air in much the same way as it occurs in 
gas-fluidized static beds. This mechanism results from an 
interplay between the flow of air through the bed and the 
effective gravitational field acting on the bed. In order 
to quantify the relevant variables of our model, we pre- 
sented explicit results for the model of Kroll. Since the 
mechanism proposed is an extension of an static one, and 
since the bed cannot switch between a fluidized state and 
a solid state instantly, we expect that this results are apli- 
cable when the period of vibration is large. Under this 
condition, we can predict when a vibrated granular sys- 
tem initially becomes gas fluidized during each cycle of 
oscillation. The proposed model constitutes a promising 
quantitative framework that gives some relevant informa- 
tion about the process of cyclic fiuidization that occurs 
in a vibrated granular bed |15|| . It would be important 
to monitor experimentally the evolution of the granular 
medium, within one period of oscillation of the system. 
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